Abstract. We map a 50' x 30' area in and around the M17 molecular complex with the 
French submillimeter balloon-borne telescope PRONAOS, in order to better understand 
the thermal emission of cosmic dust and the structure of the interstellar medium. The 
PRONAOS-SPM instrument has an angular resolution of about 3', corresponding to a 
size of 2 pc at the distance of this complex, and a high sensitivity up to 0.8 MJy/sr. 
The observations are made in four wide submillimeter bands corresponding to effective 
. . . wavelengths of 200 fim, 260 fim, 360 fim and 580 /im. Using an improved map-making 
method for PRONAOS data, we map the M17 complex and faint condensations near the 
dense warm core. We derive maps of both the dust temperature and the spectral index, 
P5 which vary over a wide range, from about 10 K to 100 K for the temperature and from 
1-^ about 1 to 2.5 for the spectral index. We show that these parameters are anticorrelated, 
the cold areas (10-20 K) having a spectral index around 2, whereas the warm areas have 
,; a spectral index between 1 and 1.5. We discuss possible causes of this effect, and we 
^ ' propose an explanation involving intrinsic variations of the grain properties. Indeed, to 
match the observed spectra with two dust components having a spectral index equal to 
2 leads to very large and unlikely amounts of cold dust. We also give estimates of the 
column densities and masses of the studied clumps. Three cold clumps (14-17 K) could 
■ be gravitationally unstable. 

o 

> 



A&A manuscript no. 

(will be inserted by hand later) 

Your thesaurus codes are: 

missing; you have not inserted them 

Submillimeter dust emission of the M17 
complex measured with PRONAOS 

X. Dupac\ M. Giard\ J.-P. Bernard^-^, N. Boudet\ J.-M. Lamarre^, C. 
Meny^ F. Pajot^, E. Pointecouteau\ I. Ristorcelli\ G. Serra\ B. Stepnik^, 
J.-P. Torre^ 

^ Centre d'Etude Spatiale des Rayonnements 

9 av. du colonel Roche, BP4346, F-31028 Toulouse cedex 4, France 
^ Institut d'Astrophysique Spatiale 

Campus d'Orsay Bat. 121 

15, rue Clemenceau, F-91405 Orsay cedex, France 
^ LERMA, Observatoire de Paris 

61, Avenue de I'Observatoire, F-75014 Paris, France 
* Service d'Aeronomie du CNRS 

BP3, F-91371 Verrieres-le-Buisson cedex, France 

Received 15 January 2002 /Accepted 12 June 2002 



Key words: dust — infrared: ISM: continuum — ISM: clouds — ISM: individual (M17) 



1. Introduction 

The submillimeter domain is particularly suited to characterizing dust properties in the 
interstellar medium. Dust emission in this spectral range is mainly due to big grains 



at thermal equilibrium (see, e.g., Desert et al. 199C), whose emission is usually mod 



elled by the modified blackbody law, i.e. by the temperature and the spectral index 
of the dust. The temperature of a molecular cloud is a key parameter which controls 
(with others) the structure and evolution of the clumps, and therefore, star formation. 
Thus spectral imaging of molecular clouds can provide useful information about their 
structure and evolution, especially if the dust emission parameters can be properly de- 
rived on top of submillimeter intensities. Mapping of star-forming molecular clouds, as 
well as other dusty regions, has been performed by the PRONAOS balloon-borne exper- 



iment (PROgramme NAtional d'Observations Submillimetriques, see Serra et al. 2001 



or Ristorcelli et al. 1998). We present in this article the maps and analysis made from 



PRONAOS observations of the M17 star-forming complex. 
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The Messier 17 Nebula (also called Omega, Horseshoe or Swan Nebula) is an ion- 
ized region associated with a giant star-forming molecular cloud located at about 2200 



parsecs from us (Chini et al. 1980) in the constellation of Sagittarius. This nebula has 
the largest known ionization rate in the Galaxy for star forming regions (see for ex- 



ample Glushkov 1998). This large ionization rate is mainly due to the excitation of gas 
by young 0-type stars. The M17 molecular cloud has been mapped in carbon monox- 
ide emission by Lada (1976), showing the two most intense condensations in this cloud, 
usually called M17 North (N) and M17 Southwest (SW). This cloud is part of a giant 
molecular complex extending 170 pc to the southwest, along the Sagittarius spiral arm 



(CO emission shown by Elmegreen et al. 1979). The interaction of this giant molecular 
cloud with the H II regions is particularly visible in the M17 region, where an expanding 
shock front interacts with the gas clouds, and is thought to have fragmented the original 



molecular cloud (Rainey et al. 1987). The M17 SW cloud is the best studied region of 
the area, especially for the photon-dominated region (PDR) near the boundary of the 
H II region (to the northeast). The first far infrared observations of M17 were made by 
Low & Aumann (1970) and Harper & Low (1971). M17 SW was mapped in the mid 
and far-infrared by Harper et al. (1976) and Gatley et al. (1979). Wilson et al. (1979) 
mapped the whole M17 cloud {i.e. including M17 North) at 69 /im. The IRAS satellite 



(http://www.ipac.caltech.edu/ipac/iras/iras.html) provided infrared maps of the whole 
sky, which showed the distribution of the dust in the M17 complex, however without 
giving much information on the faintest and coldest regions of the complex. The recent 
JCMT measurement of Wilson et al. (1999) provided precise CO maps of the M17 molec- 
ular cloud, as well as the measurement of Sekimoto et al. (1999) on a larger area. In this 
context, submillimeter mapping can provide useful information about the dust proper- 
ties, especially in cold regions, and give independent estimations of the ISM masses in 
this region. For these reasons, we have observed a large region including the M17 cloud 
and fainter areas, with the multi-band photometric instrument (SPM) of the PRONAOS 
balloon-borne experiment. The maps are 50' by 30' (about 30 pc x 20 pc at the distance 
of M17) at angular resolutions between 2' and 3.5'. We present in Section 2 how the 
observations were made and the way we processed the data to make maps. Section 3 
presents the resulting maps, and Section 4 presents a detailed analysis of these maps. 



2. Observations and data processing 

PRONAOS (PROgramme NAtional d'Observations Submillimetriques) is a French 
balloon-borne submillimeter experiment. The 2 m telescope is described in detail in Buis- 
son & Duran (1990). The focal plane instrument SPM (Systeme Photomctrique Multi- 
bande, see [Lamarre et al. 1994 ) is composed of a wobbling mirror, providing a beam 
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switching on the sky with an ampHtude of about 6' at 19.5 Hz, and four bolometers 
cooled at 0.3 K. They measure the submillimeter flux in the spectral ranges 180-240 /im, 
240-340 //TO, 340-540 /im and 540-1200 fim, with sensitivity to low brightness gradients 
of about 4 MJy/sr in band 1 and 0.8 MJy/sr in band 4. The effective wavelengths are 
200, 260, 360 and 580 fim, and the angular resolutions arc 2' in bands 1 and 2, 2.5' in 
band 3 and 3.5' in band 4. Details about the instrument can be found in Ristorcelli et 
al. (1998). The data which we present here were obtained during the second flight of 
PRONAOS in September 1996, at Fort Sumner, New Mexico. 

The usual reconstruction method from chopped PRONAOS data was EKH-like 



( [Emerson et al. 19791) with a scan by scan filtering in Fourier space (see pales et al. 199l| ). 
Dupac et al. (2001) developed another method, based on direct linear inversion on the 
whole map, using a Wiener matrix. This map-making process takes into account the beam 
sizes and profiles, the beam switching and the signal and noise properties to construct 
an optimal map. 

For these M17 data, we have used a slightly improved method, in which we consider 
some noise not independent of the sky signal. To assume a perfect independence between 
noise and signal is justified for Cosmic Microwave Background maps, in which the inten- 



sity contrasts are low (see map- making methods in Dupac & Giard 2002). However, in 
the case of timelines having much contrast such as the PRONAOS ones, the instrument 
response can hardly be considered as perfect. Thus, we allowed the noise covariance ma- 
trix to take into account noise proportionally correlated to the signal. The fraction of the 
noise that appears to be correlated to the signal is determined by iterations and tests 
of the reconstruction method. We have compared PRONAOS M17 maps made with the 
method of Dupac et al. (2001), i.e. assuming that the noise is independent of the signal, 
and the ones made with the new method. In the case of the independent-noise method, 
it appears that the found noise level constrains the lowest signal that can be recovered. 
The problem is that if the noise is stronger in intense areas, due to its correlation with 
the signal, then the independent-noise method removes signal features lower than this 
noise level. With the new method, the problem is solved by allowing the noise not to be 
the same in faint areas and in intense areas. We found that the maps made with the new 
method are significantly better as they reconstruct more faint areas. The characteristics 
of this new method are otherwise common to the one described in Dupac et al. (2001). 



3. Results 

We present in Fig. |l| the images obtained in the two extreme photometric bands of 
PRONAOS-SPM. The angular resolution is 2' in the 200 jim map, which corresponds 
to 1.3 parsec at the distance of M17, and 3.5' in the 580 fim map (2.2 pc). Due to 
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Table 1. Equatorial coordinates and fluxes (Jy) integrated over a 3.5' area, around the 
intensity peaks of the identified regions. The absolute errors on the fluxes are 8 % (1 tr). 
These values take into account the fine color corrections due to the large bandwidths of 
the SPM instrument. 
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(h,min,sec) 


5l950 

r/) 


200 fim 


260 fim 


F.(Jy) 
360 fj,m 


F.(Jy) 
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M17 SW 


18 17 35 


-16 15 


28000 


15000 
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M17 N 


18 17 45 
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Cloud A 


18 16 22 
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25 


Cloud B 


18 16 16 
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36 


Cloud C 


18 16 29 


-16 28 
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98 


31 


Cloud D 


18 15 34 


-16 11 


240 


210 


140 


44 



the calibration uncertainty, the fiux accuracy is 5 % (1 cr) relative between bands (8 % 
absolute). We also present in Table |l|the fluxes of the main identified regions, integrated 
over a 3.5' beam. 

The maps in Fig. ^ exhibit the very high intensity contrasts that exist in this kind of 
giant molecular complex. The Ml? Southwest (SW) and North (N) areas show up as the 
most intense dust emission on our maps. The Ml 7 SW cloud reaches a peak intensity of 
46000 MJy/sr in band 1 (200 /im with a 2' angular resolution) at the MlTa source (see 
Lada 197^ , Wilson et al. 1979 and Gatley et al. 1979 ). A small intense condensation is 
visible to the northeast of the most intense area. It shows up clearly in Fig. |^ with the 
second highest contour in the 200 /im image, and has a peak intensity of 9600 MJy/sr at 
200 /im. This condensation corresponds to the M17b source. The M17 North condensation 
(also called M17c) reaches a peak intensity of 7800 MJy/sr at 200 jim. The M17d source 



( Wilson et al. 1979 ) is clearly visible to the east from M17b, and reaches a 200 jim peak 
intensity of 2000 MJy/sr . 

We have also mapped the area to the west of the M17 complex, in which we can 
identify four condensations with weak intensities. Cloud A (see Fig. |^ for the names) is 
situated west of M17 SW and shows up as an extended condensation reaching a 200 /im 
peak intensity of 630 MJy/sr . Cloud B is south of Cloud A and is linked to it by a bridge 
clearly visible on the maps. Cloud B reaches a 200 jim peak intensity of 440 MJy/sr . 
Another extended condensation, that we call Cloud C, is visible south of Cloud B. It has 
a peak intensity of 470 MJy/sr . Cloud D is a condensation situated to the northwest of 
the maps, linked to Cloud B through a fllament, inside which a small (2') condensation 
is visible. Cloud D has a 200 /im peak intensity of 470 MJy/sr . 
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It is worth mentioning that these weak-intensity clouds are poorly visible in the IRAS 
100 fim survey. As we can see in Fig. |l], the PRONAOS 580 /im map looks smoother than 
the 200 /iTO one, exhibiting less intensity contrast between regions. As we shall discuss 
later, this is due to a trend of the intense areas to have warmer dust than the faint ones. 
This induces a stronger relative intensity of the weak-intensity areas in long-wavelength 
bands, therefore reducing the contrast between the regions in the long-wavelength maps. 



4. Analysis 

4-.1. Dust temperatures and spectral indices 
4.1.1. Derivation 

We assume that the emission of the grains is characterized by the modified blackbody 
law: 

I,fu{X,T,(3)^C B,{X,T)X~I' (1) 

where A is the wavelength, C a constant, T the temperature of the grains, (3 the 
spectral index and the Planck function. The fit procedure is described in detail in 
Dupac et al. (2001). In this work, we included iteratively the color correction process 
due to the SPM bandwidths. Assuming optically thin emission, we can derive the optical 
depth = ^ = C X-f^ . The relative error on r,, is + -^AT. 

We show the temperature and spectral index maps derived this way in Fig. |^. These 
have been made without assuming anything about the temperature or the spectral index. 
We used the IRAS 100 /im data for the determination of the temperature and the spectral 
index of a large part of the map (for this we assumed an intercalibration error between 
IRAS and PRONAOS data of 25 %), but neither for M17 SW (saturated) nor for the 
faint clouds A, B, C, D to the west, for which the IRAS 100 /im data seem too noisy. In 
the areas where the fit was difficult because of the high temperatures estimated (> 70 
K), we used the IRAS 60 fim data too. This is justified by considering that the thermal 
continuum of big grains dominates the 60 /zm emission for such high temperatures, 
whereas it is likely not to be the case for lower temperatures for which the very small 
grains might dominate the 60 /im emission. With this process, the fit is performed in 
each pixel where the signal to noise ratio is reasonable. However, some observed areas 
are too noisy to obtain an estimation of the temperature and the spectral index. These 
areas are displayed in white on the maps in Fig. ^. We present spectra of some observed 
Ml 7 regions in Fig. |. 
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Table 2. Temperature, spectral index and optical depth of the intensity peaks of the 
identified regions. The error bars are given for the 68 % confidence interval. 





T (K) 




r,. . 10^ 

200 ^^m 


. 10^ 

580 /ijn 


M17 SW 


29 ± 8 


1.7 ± 0.3 


80 ± 60 


12 ± 3 


M17 N 


28 ± 3 


1.6 ± 0.2 


15 ± 6 


2.6 ± 0.4 


Cloud A 


17 ± 3 


2.3 ± 0.3 


6 ± 4 


0.5 ± 0.1 


Cloud B 


17 ± 3 


1.7 ± 0.3 


5 ± 3 


0.7 ± 0.2 


Cloud C 


26 ± 4 


1.2 ± 0.2 


1 ± 0.5 


0.3 ± 0.07 


Cloud D 


14 ± 2 


1.9 ± 0.3 


9 ± 6 


1.3 ± 0.3 



4.1.2. Variations of the temperature and the spectral index 

The maps in Fig. exhibit interesting features: first, it is evident comparing the maps 
that an anticorrelation exists between the temperature and the spectral index. We shall 
quantify this later. In the major part of the maps, the temperature ranges from 10 K to 
80 K, while the spectral index exhibits also large spatial variations from 1 to more than 

2.5. Temperatures above 80 K (around 100 K or higher) can be found in a few pixels near 
the ionization front to the east. Large scale features on these maps are clearly visible, 
showing how the temperature increases from the west to the east. The weak intensity 
areas to the west are significantly colder than the intense M17 area. We present in Table ^ 
the temperature, spectral index and optical depth for the intensity peaks of the identified 
regions (Fig. |l]). The M17 SW region shows a temperature of 29 K (± 8) at the intensity 
peak. For this area, we derive a spectral index of 1.7 (± 0.3). However, since the IRAS 
data are saturated in this area, the fit is quite uncertain and the temperature could be 
higher. The M17 N condensation has a temperature of 28 K, and a spectral index of 

1.6. The temperature in the M17 complex is about 30-50 K outside the M17 SW and 
M17 N condensations, and the spectral index varies between 1 and 1.5. On the edges of 
the complex, one can see the warm areas near the ionized regions (east and north-west), 
where the temperature can be higher than 80 K, and the spectral index between 0.7 and 
1.1. The clouds A, B, and D in the western part of the maps show low temperatures 
(14-17 K) and high spectral indices (1.7-2.3). Some areas in these clouds even exhibit 
very low temperatures down to 10 K. 

4.1.3. Anticorrelation between the temperature and the spectral index 

As we can see in Fig. |^, it seems that an anticorrelation exists between the temperature 
and the spectral index. This effect was shown for the first time in observations by Dupac 
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et al. (2001) in Orion. If we take the five intensity peaks presented in Table ||, the 
correlation coefficient is -0.58, which means a significant anticorrelation. However, it is 
not highly statistically significant, given the small number of points. Therefore, we have 
analyzed the correlation between the temperature and the spectral index in the whole 
map. We present in Fig. ^ a plot of the distribution of the pixels in Fig. ^ in the (T,/3) 
space. 

One can see the banana shape of this distribution in Fig. ^ which clearly shows the 
anticorrelation between both parameters. In particular, very few pixels can be found 
in the data with T > 40 K and (3 > 1.5, and no pixel can be found with T < 20 
K and /3 < 1.6. If we consider all the (T,/3) pairs with T between 10 K and 80 K, 
and f3 between 1 and 2.5, then the correlation coefficient is -0.56. Among these 1079 
pixels, 1071 have both temperature and spectral index with relative errors less than 50 
%, and 801 with relative errors less than 20 %. The correlation coefficient computed 
from these 801 well fitted pixels is -0.64. If we consider only the pixels with relative 
errors less than 10 % (87 pixels), then the correlation is -0.74. This increasing of the 
anticorrelation with the restriction on the errors is evidence of the reality of this effect. 
Also, pixels with high spectral indices have narrow error bars on the spectral index: 
there are 70 pixels for 2.3 < /3 < 3, with an average relative error on /3 of 12 % - 
the relative errors are quite uniform for these points. Pixels with low spectral indices 
also show quite narrow error bars: there are 482 pixels for 0.7 < f3 < 1.3, with an 
average relative error on (3 of 16 %. Therefore, we can be confident of this anticorrelation 
effect. Nevertheless, to better understand it, we performed simulations of PRONAOS data 
from the same number (801) of random (T,/3) pairs in the same (T,/3) range, and fitted 
them the same way as real data. We repeated this processing with several distributions 
of T-/3 pairs without intrinsic correlation, with and without 100 /im simulated data, 
and obtained correlation coefficients on the fitted parameters usually between 0.0 and 
-0.2, and never above (in absolute) -0.3. Therefore, the fit procedure itself involves some 
degree of anticorrelation between the derived temperature and spectral index, as Dupac 
et al. (2001) had already investigated. There is indeed some degeneracy between both 
parameters, but this correlation coefficient which we find on simulations is clearly not 
enough to explain the correlation found on the data (-0.64). Thus we have to conclude 
that this anticorrelation effect has to be either an intrinsic physical property of the grains, 
or at least a property of the observed dust, when integrated on the beam and the ISM 
column. The same effect has been shown in Orion ( pupae et al. 200"l ). 

To interpret this anticorrelation, we can investigate what the meaning is of such 
spectral index variations. Mixtures of dust components with different temperatures can 
be put forward: this clearly cannot explain high indices, but could explain low indices 
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Table 3. Summary of the investigations made concerning the two-dust component as- 
sumption. From left to right: the original spectrum temperature, the original spectral 
index, the temperature of the cold component, which has a spectral index (3=2, and the 
mass ratio derived between the cold and warm components. 



T (K) 


/3 


Tcoid (K) 


Mass ratio 


50 


1 


10 


100 


70 


1 


10 


200 


60 


1 


20 


20 


60 


1 


40 


no fit 


50 


1.5 


10 


50 


50 


1.5 


20 


5 


50 


1.5 


30 


3, bad fit 


40 


1.7 


10 


20 



around 1. Indeed, warm dust can be associated in the same line of sight with cold dust 
components which could enhance the high-wavelength intensity, and therefore decrease 
the spectral index from 2 to around 1. However, this would imply a large amount of 



cold dust (see the investigations of Ristorcelli et al. 1998 and Dupac et al. 2001). In this 
article, we have performed other studies. We simulate spectra with low spectral indices 
and various warm temperatures. We then try to fit these spectra with two dust compo- 
nents having both a "standard" spectral index equal to 2, the warm component having 
the same 100 /^m emission and the same temperature as the single component spectrum. 
Thus we have equality between both optical depths (C A"'') at 100 jim, and therefore. 



given the simple column density model that we present in Section 4.2, the warm dust 
component has the same column density as the single dust component. In this way, we 
fit the short-wavelengths part of the original spectrum by the warm component, but the 
discrepancy between the original spectral index and the standard one (/3 = 2) induces 
a lack of emissivity of the warm component compared to the original spectrum at long 
wavelengths. Therefore, we introduce a colder dust component with the same standard 
spectral index /3=2, which we add to the warm component to fit the spectrum. Then we 
can investigate what amount of cold dust is necessary to explain the observed low spectral 
indices, by computing the column densities from our simple model (see Section [I.2| and 
Dupac et al. 2001). Indeed, with two dust components with the same standard spectral 
index, the column density ratio between both is simply the ratio of the parameters C. 

We present in Table |^ a summary of these investigations. It appears that low indices 
(around 1) would imply very large amounts of cold dust in the line of sight. For example, 
to fit this way a spectrum with T = 50 K, /3 = 1, with a 10 K component, a cold dust 
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column density 100 times larger than the warm dust column density is necessary. We 
have said that the column density of the warm component is the same as the single 
dust-mcthod one. As seen in Table |4| the column density estimates, assuming one dust 
component, seem to have the right order of magnitude, compared for example to the 
estimates from the CO. It is extremely unlikely that there are such huge amounts of cold 
dust (much larger than the warm dust and the CO estimates) in each line of sight where 
we observe low indices. If one uses intermediate temperatures for the cold component, 
then it is no longer possible to match the spectrum at long wavelengths, even for 
We show by these investigations that even for intermediate spectral indices, it is unlikely 
that they be explained only by temperature mixtures with the same standard spectral 
index. 

Moreover, the shape of the pixel distribution in the (T, /3) space is very difficult 
to explain with only the argument of temperature mixtures. It would suggest that the 
warmer the (warm) dust component, the more massive the cold component, because the 
spectral index is further reduced (to 1) at high temperatures. This can certainly not be 
a general rule, and this shows that the temperature mixture assumption is insufficient 
to explain the variations found on the spectral index and its inverse dependence on the 
temperature. 



Laboratory experiments (see Agladze et al. 1996 and Mennella et al. 1998 ) showed 
this anticorrelation effect on grains for temperatures down to 10 K. Agladze et al. (1996) 
measured absorption spectra of crystalline and amorphous grains between 0.7 and 2.9 
mm wavelength. They deduced an anticorrelation between the power-law index [3 and the 
temperature in the temperature range 10-25 K, and attributed it to two-level tunnelling 
processes. The measures of Agladze et al. (1996) are insufficient to justify our observation 
in the submillimeter spectral range, because absorption can be very different in the 
millimeter range. Mennella et al. (1998) measured the absorption coefficient of cosmic 
dust grain analogues, crystalline and amorphous, between 20 /zm and 2 mm wavelength, 
in the temperature range 24-295 K. They deduced an anticorrelation between T and /3, 
and attributed it to two-phonon difference processes. In our observations, we observe 
this effect down to about 10 K (Fig. thus we would need laboratory results on these 
low temperatures in the submillimeter range to fully understand the observations. Other 
causes can make the spectral index vary, such as the composition and size of the grains. 

We can conclude from this temperature-spectral index analysis that the anticorrela- 
tion discovered by Dupac et al. (2001) in Orion is shown in M17 too, another high-mass 
star-forming complex of our Galaxy. Our investigations make us rather believe in a fun- 
damental property of the grains to explain this effect. 
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Table 4. Column densities estimated from the PRONAOS+Desert method, and 
from the ^^CO data of Wilson et al. (1999), masses from PRONAOS+Desert, 
PRONAOS+Osscnkopf, Jeans masses. The error bars on the PRONAOS+Desert col- 
umn densities are given for the 68 % confidence interval. 





N/f PRONAOS + Dgs. 




Mass PRONAOS + Doa. 


Mass PRONAOS + Osa. 


Jeans mass 




lO^^Hcm-^ 






M© 


Mo 


M17 SW 


4100 ±3600 


1200 


16000 






M17 N 


730 ±340 


460 


4600 






Cloud A 


450 ±380 




2000 


730 


70 


Cloud B 


250 ±210 




2400 


870 


100 


Cloud C 


38 ±25 




350 




150 


Cloud D 


630 ±530 




2100 


760 


100 



4-2. Column densities 

We are interested in estimating the column densities and masses of the studied region. 
For this we follow the simple model described in Dupac et al. (2001), which uses the dust 
100 /im opacity from Desert et al. (1990). We derive column densities that we present in 
Table ^. To obtain an independent estimation of the column density, we use observations 
of the rotational transition J = 1 — of ^'^CO, made by Wilson et al. (1999). Their 
map covers the M17 cloud but not the fainter area to the west. Assuming local thermal 
equilibrium and optically thin emission, the ^^CO column density can be expressed as: 
^-co = 2.6 10-^^^^ (2) 
(see e.g. Rohlfs et al. 2000) Then, assuming a li^Q number ratio of 4.6 10^ 



(Rohlfs et al. 2000), a mean molecular weight of 2.36 m// (see e.g. Elmegreen et al. 1979), 
and a temperature of the gas of 30 K, we derive a column density to ^^CO line ratio of 
17 10^° protons cm^"^ [K .km / s)^^ . Using this coefficient, we derive column densities that 
we present in Table ^, in order to compare to the PRONAOS + Desert model estimation. 

The agreement is not very good in M17 SW between PRONAOS + Desert model and 
^•^CO, though it is consistent with the error bar. It is not unlikely that this is due to 
optical thickness of the -'^'^CO 1-0 emission. The M17 N region shows a better agreement. 
If we trust our PRONAOS + Desert model estimation (which seems quite speculative in 
view of the error bars), we can constrain the column density to ^^CO line ratio in M17, 
which we find to range roughly from 15 10^" to 100 10^" following the regions of M17 
(in units of protons cm^'^ [K .km / s)^^) . Since the ratio corresponding to optically thin 
^^CO emission is 17 lO^*', and that the most intense areas have the highest coefficient. 
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it is possible to conclude that the disagreement between PRONAOS and ^^CO column 
density estimations of M17 SW may be due (at least in a part) to optical thickness of the 
^^CO 1-0 line. Also, uncertainties in the ri^Q ratio may arise for high column densities. 
Moreover, the excitation temperature can be higher than 30 K in M17 SW, although 
we measure it to be around 30 K. If we take an excitation temperature of 60 K, then 
the column density estimated from the ^■^CO emission is 2300 lO^*' H cm^^ for M17 SW 
(at the intensity peak). If now we look at the sources of uncertainties in the PRONAOS 
+ Desert model estimation - apart from the error bars being large, we can point out 
that the 100 iim opacity we used might overestimate the column density, especially in 
cold clouds where could take place some special effects like formation of molecular ice 
mantles on the grains and coagulation of grains. For instance, these processes are taken 
into account by the protostellar core model of Ossenkopf & Henning (1994), for which the 
100 opacity is 1 cm? / g. By making the same analysis as with the model of Desert et 
al. (1990), using the Ossenkopf & Henning opacity, we derive column densities 2.8 times 
lower. This correction could be true for cold clouds like A, B and D. 

Finally, we computed masses of regions in this giant molecular complex. For this, 
we integrated the column density found from our submillimeter measurement over 
the area of the clouds that we observe in the maps, assuming a distance of 2200 pc 
( |Chini et al. 1980 ). By this way, we derive a total ISM mass of the M17 complex (with- 
out Clouds A, B, C, D to the west) of 31000 M©, in which M17 SW weights 16000 M© 
and M17 N 4600 M©. We derive also mass estimates of the clouds to the west, which we 
present in Table ^. For the cold clouds A, B and D, it is likely that the masses are smaller 
than that: if we trust the Ossenkopf & Henning (1994) model, then we derive masses of 
about 800 M0 for each. For the clouds separated from the molecular complex (A, B, 
C, D), we present also in Table ^ the Jeans masses, as an indicator of the gravitational 



stability of these condensations (see [Dupac et al. 200 1| for details). There are of course 
large uncertainties in both the clump's measured mass and the non-thermal sources of 
internal energy, such as magnetic field and turbulence, which can help to balance the 
gravitational energy. Nevertheless, we observe that Clouds A, B and D are possibly grav- 
itationally unstable, especially Cloud A, which might be compressed by the external 
pressure due to the ionized area to the north. 



5. Conclusion 

We observed a large (50' x 30', w 30 pc x 20 pc) area in and around the M17 molecular 
complex. We showed extended faint clumps (A, B, C, D) in the outskirts of the active 
star-forming area. Our study shows a large distribution of temperatures and spectral in- 
dices: the temperature varies roughly from 10 K to 100 K, and the spectral index from 1 
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to 2.5. The statistical analysis of the temperature and spectral index spatial distribution 
shows the anticorrelation between these two parameters. Indeed, we observe cold dust 
(10-20 K) with high indices (around 2), and warm dust (> 20 K) with low indices (around 
1 - 1.5), but we could not significantly find any cold place with low indices nor warm 
place with high indices. The investigations which we made to match the observed spec- 
tra, with two dust components having a standard spectral index, imply very large and 
unlikely amounts of cold dust, therefore we rather support a fundamental explanation 
for this effect. Laboratory measurements showed this anticorrelation in the submillimeter 
domain on grains for temperatures down to 25 K, but we would need laboratory results 
on temperatures down to 10 K in the submillimeter range to fully understand our obser- 



vations. This anticorrelation effect was also shown in Orion (Dupac et al. 2001), and in 
other PRONAOS observations that are still being analyzed. 

We estimated the column densities and masses of the observed regions by simply 
modelling the thermal emission of the grains from Desert et al. (1990) or Ossenkopf & 
Henning (1994). We derived a total mass of the M17 cloud of 31000 Mq. Of course, 
the error bars are large and difficult to estimate properly, but there is a clear trend for 
three cold clouds (A, B, D) to be gravitationally unstable, and therefore to be pre-stellar 
candidates. 

These observations could be sustained by sensitive continuum observations in the 
millimeter domain, especially for the faint clouds studied to the west of the M17 complex, 
in order to better constrain the spectral index measurement of the cold dust. Also, 
infrared observations at higher resolution could be useful to study the structure of these 
cold clouds. 
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Fig. 1. PRONAOS intensity maps in band 1 (200 /im, up) and band 4 (580 fim, bottom). 
The angular resolution is 2' for the 200 ^im map and 3.5' for the 580 /im map. The color 
scale is logarithmic, and displays the positive reconstructed intensity until -1 in log, then 
the deeper blue and purple colors display the negative noise features. The noise level 
is about 4 MJy/sr rms in band 1 and 0.8 MJy/sr in band 4. Due to the calibration 
uncertainty, the intensity accuracy is 5 % (1 ct) relative between bands (8 % absolute). 
The non-observed areas are displayed in white. 
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Fig. 2. Spectra with 3 a error bars. The PRONAOS error bars are the intcrcahbration 
errors. The 100 fxm point is from the IRAS survey. The drawn hnes are the result of the 
fits (modified black body). 
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Fig. 3. Maps of the dust temperature in Kelvin (top) and the spectral index (bottom), 
made by fitting the modified blackbody to the PRONAOS spectra, with the IRAS 100 
fim and 60 fim data for some regions. The angular resolution is 3.5' in each map. The 
external contour shows the limit of the observed area. The white regions inside this 
contour are too noisy to obtain an estimation of the temperature and the spectral index. 



18 



Dupac et al: Dust emission in M17 



3.0 



2.5 



2.0 



1 .5 



1.0 



0.5 I , , , I , , , I , , , I i , ^ J , , , 

20 40 60 80 100 

Temperoture (K) 



Fig. 4. Distribution in the (T,/3) space of the (T,/3) pairs derived with both errors on 
and /3 lower than 20 %. 



